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Abstract 

Small-world (SW) networks have been identified in many different 
fields. Topological coefficients like the clustering coefficient and the 
characteristic path length have been used in the past for a qualitative 
characterization of these networks. Here a dynamical approach is used 
to characterize the small-world phenomenon. Using the /3— model, a 
coupled map dynamical system is defined on the network. Entrance to 
and exit from the SW phase are related to the behavior of the ergodic 
invariants of the dynamics. 

Keywords: Small World, Ergodic invariants. Coupled maps. Phase tran- 
sitions 



1 Introduction 

Networks are prevalent in all domains of life and science. Social, economic 
and political networks are the backbone of human society. The internet is a 
network. The metabolic processes of living beings form a network with the 
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substrates as nodes, which are hnked whenever they participate in the same 
biochemical reaction. Protein-protein and gene expression and regulation are 
biological networks, etc. 

Regular lattices and random graphs [|l| have for a long time been stud- 
ied. More recently 0-0 small- world networks became the object of growing 
attention and were identified in many different fields. They seem to be the 
underling structure for some important phenomena like the rapid spread of 
diseases, social networks, cooperative behavior between competing agents, 
problem solving organization and communication networks. 

Topologically, small-world (SW) networks are identified by the values of 
two statistical properties: 

• the clustering coefficient (CC) that measures the average probability 
for two agents, having a common neighbor, to be themselves connected 
and 

• the characteristic path length (PL), this being the average length of the 
shortest path connecting each pair of agents. 

Regular lattices have long path lengths and high clustering, whereas ran- 
dom graphs have short path lengths but low clustering. SW networks exhibit 
short PL's and, at the same time, high C(7s. 

In many model networks, the simultaneous occurrence of high CC and 
low PL is observed over an interval between order and randomness, which 
is called the SW phase. However, this phenomenon can only be defined as 
a phase, in the statistical mechanics sense, if order parameters are found to 
characterize the regular-to-SW and the SW-to-random phase transitions. 

Further information on the SW phenomenon has been obtained in the past 
from the study of several quantities. Farkas et al. 0] studied the spectral 
density of the adjacency matrix, with increasing randomness, concluding 
that, in spite of the blurring of singularities, a consistently high value of the 
third moment implies the existence of a large number of triangles in the SW 
network. Monasson |^ on the other hand studied the spectral properties of 
the Laplacian operator, that characterizes the time evolution of a diffusive 
field, and localization properties on the graph. 

In this paper a dynamical systems approach is used to characterize the 
small- world phenomenon. Using the /3— model ||^, we study a coupled map 
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system on the network, with interactions defined by the network connections. 
The SW phase is related to the behavior of the ergodic invariants of the 
dynamics. Entrance to the SW phase is related to the Lyapunov spectrum 
and exit from the SW phase corresponds to the region where "entropy" and 
"conditional exponents entropy" split apart. 



2 The dynamical model 



Consider a /3— family of models, each one with agents on a circle and pe- 
riodic boundary conditions. For /3 = 0, each agent in the model is connected 
to its 2v nearest neighbors. For jS ^ 0, the network structure is obtained 
by looking at each one of the connections of the /3 = structure and, with 
probability f3, replacing this connection by a new random one. 

On each one of the networks, a dynamical system is defined, with a 
map at each node and convex-coupling interactions defined by the network 
connections 

N 

Xi{t + l) = Y.W,,f{x^{t)) (1) 



where 



I-^c if z = j 
~ { ^ if ^ 7^ J and i is connected to j (2) 

otherwise 

riv (i) is the number of agents connected to i and c is a control parameter. 
For the agents dynamics we choose 

/ (x) = ax mod. 1 (3) 

Typically a = 2. 

For the (3 = network, each agent has exactly 2v neighbors and the 
Lyapunov exponents are 

Ao (fc) = log |a - c + ^ ^ cos (jOk) ] \ (4) 

with 9k = ^jf , k = 0, ■ ■ ■ , N — I. In the N ^ oo limit, the Lyapunov 
spectrum is a continuous smooth function, as illustrated in the upper plot of 
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Figure 1: Lyapunov spectrum for /3 = and for a typical network at /? = 0.2 
{N = 800, 2v = 6) 

Fig. 1. For reasons to be discussed below we always choose c in such a way 
that, ior P — 0, the lowest Lyapunov exponent is zero. 

As /? increases, the matrices of the tangent map cease to be regularly 
organized, the Lyapunov spectrum develops gaps and some of the exponents 
become negative. This is illustrated in the lower plot of Fig. 1 for N = 800 
and 2v — 6. 

It is also the appearance of random long range connections that is re- 
sponsible for the reduction of the path length in SW networks. Therefore 
it is natural to consider the modifications in the Lyapunov spectrum as the 
dynamical signature of the onset of the SW phase. Of particular dynamical 
significance is the shift of part of the spectrum towards negative values. That 
is, the randomness arising from the rewiring leads to an effective reduction 
of the dynamical degrees of freedom. We define Dp 

Dp^-^X^ (5) 

Ai<0 

to quantify this effect. To characterize the modifications of the Lyapunov 
spectrum, another possibility would be to measure the singular part of the 
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spectrum associated to the gaps. However the natural intervals in the spec- 
trum associated to finite make this measurement less reliable. 

In the upper plot of Fig. 2 we show the average values of D/j taken over 
100 different samples for each (3 (with = 800 and 6 as the average degree 
of the network). A good fit to all the data shown in the lower plot of Fig. 2 is 

Dp = cN {(3 - p,,y' (6) 

with < 10"^ and r]i = 1.01 ± 0.06. 
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Figure 2: The parameter ("SW order parameter") averaged over 100 
sample networks (A^ = 100, ■ ■ ■ , 800, 2v = 6) 

In practice it is only after (3 ~ 10~^ that small-world effects (and 
values) become appreciable. Nevertheless, the fact that the data is consistent 
with = implies that, using j^D/^ as an order parameter for the small- 
world phase, this phase starts at /3 = 0"*", the regular phase being only the 
isolated point (3 = 0. 

To characterize the exit from the SW phase, we use the notion of con- 
ditional Lyapunov exponents . They were introduced by Pecora and Carroll 
in their study of synchronization of chaotic systems^]. Like the Lyapunov 
exponents, the conditional exponents are well defined ergodic invariants H]. 
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The idea is that the conditions that in Oseledec's theorem insure the exis- 
tence of the Lyapunov exponents also estabhsh the existence of characteristic 
exponents formed by subblocks of the tangent map matrix. Here, for each 
agent i, we consider a subblock of dimension di x di formed by himself and 
those that are connected to it. The positive conditional exponents A^(j') 
associated to each subblock are computed and a dimension-weighed sum is 
performed over all subblocks. This gives a version of what elsewhere [^] 
has been called a conditional exponents entropy. 



N 



( 
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i=i y"-* A*>o 



(7) 



Subtracting from the sum of the positive Lyapunov exponents, h/^ 
J2xa>o ^pU)^ "W'e define the coefficient 



h* - ho 



h* 



hr- 



(8) 



which is also an ergodic invariant. 

This coefficient has the following dynamical interpretation: The Lya- 
punov exponents measure the rate of information production or, from an 
alternative point of view, they define the dynamical freedom of the system, 
in the sense that they control the amount of change that is needed today to 
have an effect on the future. In this sense the larger a Lyapunov exponent 
is, the freer the system is in that particular direction, because a very small 
change in the present state will induce a large change in the future. The con- 
ditional exponents have a similar interpretation concerning the dynamics as 
seen from the point of view of each agent and his neighborhood ||^ . However 
the actual information production rate is given by the sum of the positive 
Lyapunov exponents, not by the sum of the conditional exponents. There- 
fore, the quantity h*^ — hp is a. measure of apparent dynamical freedom (or 
apparent rate of information production). As self-organization in a system 
concerns the dynamical relation of the whole to its parts, this quantity may 
also be looked at as a measure of dynamical selforganization. 

In Fig. 3 we show the average values of Cp taken over 100 different 
samples for each (3 (with 6 as the average degree of the network and = 
100, 200, 400, 600, 800). Notice the A^-independence of Cp which follows from 
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Figure 3: The parameter ("SW exit") averaged over 100 sample networks 
{N = 100,- ■■,800, 2v = 6) 

the fact that, in Eq.@ it is defined as a ratio of two quantities with the 
same A^-dependence. For small /3 values the difference between the entropy 
and the conditional exponents entropy is a small quantity, that may be eas- 
ily computed from the network parameters. It means that each agent may 
have exact information on the global behavior from observation of his own 
neighborhood. When (3 increases the difference changes sign and becomes 
very large meaning that the neighborhood information has ceased to provide 
reliable information on the actual dynamics. This is the dynamical correlate 
of the decreasing cluster properties and allows us to define the transition at 
the divergence point /?C2 of C'/3. We find 

Pc, - 0.04 (9) 

Near the transition region 

c^^\f3-p,,r' (10) 

with r]2 — 1.14 below the transition and r]2 — 0.93 above it. 
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3 Conclusions 



1 - The ergodic invariants (Lyapunov spectrum and conditional exponents) 
provide a link between the topological properties of SW networks and the 
dynamical behavior of a coupled map system modelled on the network. In 
addition, the power laws obeyed by these invariants provide a framework to 
identify the SW phenomenon as a phase in the statistical mechanics sense. 

2 - Coupled map behavior, evolution of a diffusive field and spectrum 
of the adjacency matrix [Q supply complementary information on the SW 
phenomenon. It is therefore conceivable that quantities obtained from these 
other approaches might also be used to construct order parameters charac- 
terizing the SW phase. 
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